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Non-Fermi Liquid Behavior of U-impurity with /^-Singlet Ground State 
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It is shown by the Wilson's numerical renormalization-group approach that Non-Fermi liquid 
(NFL) behaviors arise even for U-impurity with an /^-singlet ground state, and all the puzzling 
behaviors of Ri_2,UiRu2Si2 (R=Th, Y and La, x < 0.07) can be explained on its basis. There 
exists an unstable NFL fixed point due to the balance between Kondo singlet and Crystalline- 
Electric-Field singlet ground states. For wide region in parameter space the system shows the 
NFL behavior as a transient phenomenon in temperature range accessible by experiment. Results 
for temperature dependence of the resistivity, the magnetic susceptibility and the Sommerfeld 
coefficient with and without magnetic field reproduce the behaviors observed in Thi_a;Ua;Ru2Si2. 
Less anomalous properties of Ri_iUa:Ru2Si2 (R=Y and La) can be understood along the same 
scenario provided that those compounds have different sets of parameters. 

KEYWORDS: Non-Fermi liquid, numerical renormalization-group method, Crystalline-Electric-Field singlet, 
Ri_j;Ua;Ru2Si2 (R=Th, Y and La) 



Recently, issues about the Non-Fermi liquid (NFL) be- 
haviors in dilute U arroys have attracted much attention. 
They are related to a general interest of Kondo phenom- 
ensB^ which arise in some cases due to the existence of 
plural number of strongly correlated 5f electrons per U 
ion. Especially, the behaviors in series of compounds 
Ri_;j;U;j;Ru2Si2 (R=Th, Y and La, x < 0.07) have been 
well investigated as a single-site effect, because anoma- 
lous properties of these jmaterials are well scaled by im- 
purity concentrationsBao' The latter fact makes the 
theoretical approach simple compared with the system in 



which the intersite efS 
such as Yi-icUicPda 



pt is considered to play some role, 
& Thi_.^U^Ru2Si2 shows - InT 
divergence both of the susceptibility x ^^'^ the Sommer- 
feld coefficient 7 with decreasing temperature. However, 
a big puzzle still remains as discussed shortly. 

The logarithmic dependence of these thermodynamic 
quantities itself is consistent withjihe prediction based 
on the quadrunolar Kondo modclj3^ if it is generalized 
appropriately .Bu^ The quadrupolar Kondo effect is ex- 
pected to arise for a non-Kramers doublet ground state 
of /^-configuration and to he mapped onto the S=l/2 
two-channel Kondo model JlSP The exact solution of 
S=l/2 two-channel Kondo model predicts the asymp- 
totic behavior of 7, x —InT and the resistivity p 
(X const ± Vt in the limiJ-JZLpdi^O with the residual en- 
tropy Sir -> 0) = iln2y'BB Since Ri_^U:,Ru2Si2 
(R=Th, Y and La) has the tetragonal symmetry and 
the valence of the uranium ion is expected to be mainly 
U'*+ (5/^), many theoretical studies of this problem have 
been done on the basis of the two-channel Kondo model 



assuming the magnelk 
(5/^) configuration.BBlii 
However, there exist 



round state, IFL ± 



puzzling properties about 
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Thi_a;Uj;Ru2Si2 which cannot be explained by the two- 
channel Kondo model. First of all, x shows the Curie 
law while p shows large value of the order of the uni- 
tarity hmit for 10K< T < lOOK. Namely, the suscep- 
tibility indicates the existence of the localized moment 
while the resistivity suggests the confinement of the lo- 
calized moment by the spin polarization of conduction 
electrons. This apparent inconsistency was removed by 
the extended two-channel Anderson model intmdjji^ed 
and discussed by Sakai and his collaboratorsjij'c3E3' 
who showed that there exists wide set of parameters 
which gives results consistent with the experiments of 
Thi_a;Ua;Ru2Si2. Howcvcr, in the case of Ri_j;Ua;Ru2Si2 
(R=Y and La), p is almost constant below lOK, which 
cannot be solved by such an extended two-channel An- 
derson model. Secondly, the divergence of 7 is sup- 
pressed for applied magnetic field H of the order of sev- 
eral Tesla in Thi_2;U2:Ru2Si2 Jj) In any theories based on 
the two-channel Kondo model, the entropy is expected to 

be released at finite temperature because the magnetic 

f 21 

field lifts the degeneracy of the |IlL ±) doublet ground 
state, leading enhancement of 7.EZP Thus, the theory 
based on the two-channel Kondo model with magnetic 
doublet configuration does not seem to be a final theory 
for this problem. 

The purpose of this paper is to solve this puzzle on 
the basis of a novel mechanism assuming Crystalline- 
Electric-Field (CEF) singlet ground state, in which the 
competition between the conventional Kondo singlet and 
CEF singlet gives rise to the NFL behavior in low but 
intermediate temperature regime as a transient phe- 
nomenon. A system exhibiting such competition can 
be mapped to the /^-impurity model with anisotropic 
antiferromagnetic Hund's-rule coupling and different hy- 
bridizations with conduction electrons, when CEF lev- 
els are specified by pseudo spins. It has been recog- 
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nized by recent studieflBBBBB that such mod- 
els exhibit the NFL behavior reminiscent of two-channel 
Kondo model for certain critical set of parameters. On 
the basis of Wilacm's-pumerical renormalization-group 
(NRG) method Jd' £3' we show that the NFL behav- 
iors arise, as a transient phenomenon near the unsta- 
ble fixed point, rather robustly even for parameters off 
the criticality in rather wide temperature region acces- 
sible by experiments. Series of present results are con- 
sistent with those observed in Ri_2:U2:Ru2Si2 (R=Th, 
Y and La). Especially, all the anomalous properties of 
Thi_a;UxRu2Si2 can be explained semiquantitatively for 
a certain set of parameters. 

We consider two low- lying doublet states out of 
three doublet oi j = 5/2 orbitals in the tetragonal sym- 
metry. These four (two doublet) states can be repre- 
sented by pseudo spin: 



.(2) 



.(2) 

-7 ' " 



|r6,+) 



|r6,-) = 




--)^|T,o), (0.1) 



U,0), (0.2) 
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where | t, 0) (|0, t)) represents the state that the pseudo 
spin in channel 1 (2) is up and the channel 2 (1) is 
empty, and notations are conventional otherwise. In 
terms of these four states, the relevant CEF states in 
/^-configuration are written in the j-j coupling scheme 
as follows: 
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Here, the energy levels of | tiiO) and |0, ti) have been 
omitted because of large intra-orbital Coulomb repulsion. 
One of fundamental assumptions is to adopt CEF level 
scheme shown in Fig.l. This level scheme can be repro- 
duced by the pseudo-spin Hamiltonian 

H - 
-nHund — 



-[8^82 



3^:3+] + J, SfS^, (0.9) 



where 3m denotes a pseudo-spin operator of the localized 
electrons in the orbital m. By taking such the pseudo 
spin representation, it can be understood clearly that the 
CEF singlet competes with the Kondo singlet directly 
and we can perform calculations more accurately owing 
to an appearance of new conserved quantity, pseudo spin. 
The energy separation between the singlet (r4) and the 

(2) 

doublet (Fg '±) is set to be A and that between the r4 
and the singlet (Fa) is set to be K. The couplings J± 
and Jz are related to {K, A) as J±=K and Jz=2A — K. 



It is remarked that (0.9) represents an anisotropic anti- 
ferromagnetic Hund's-rule coupling, which cannot exist 
in real spin system. 

With this preliminary, the system in question can be 
described by a conventional two-channel Anderson model 
with antiferromagnetic Hund's-rule coupling: 



Hf + Hy{ 



und ; 



with 



V Ve-ct c- 



m=l,2 r 



ka 



^ kma 



h.c). 



(0.10) 



(0.11) 



(0.12) 



— 'y ^ E[„ifma-fma + ^ ^ 2 fma-frncrfmafnJfi-^'^) 
ma ma 

It is noted that the hybridization ^ depends only 
on its channel index m{= 1, 2) and the inter-orbital 
Coulomb interaction is neglected for simplicity. Here- 
after, we take the unit of energy as (l-t-A~^)_D/2 where D 
is half the bandwidth of conduction electrons and A is a 
discretization parameter in NRG calculation. In the case 
A'=:A=G, this model reduces to two independent conven- 
tional Anderson models. Therefore, we can restrict our- 
selves to the case that each /-orbital is almost singly oc- 
cupied by taking parameters such that En — Ef2 — —0.4, 
Ui = U2 = 2.0, Vi = 0.7, and V2 = 0.34. 

The system described by (0.10) shows the conipfvtition 
between the Kondo singlet and±hp_]cK:al^ingletE§l^) as 
seen in two-impurity systems.E3EilE3 E3ES Two types of 
ground state are possible, i.e., (i) the Kondo singlet state 
characterized by total phase shift being equal to tt (i.e., 
Si = 62 = 7r/2), and (ii) the local singlet state charac- 
terized by Si ^ 62 = 0. Then, there exists a locus of un- 
stable fixed points characterized hy 61 = S2 = /4^, and 
on which the NFL behavior appears. A critical line can 
be determined by the NRG calculation, across which the 
even-odd alternation in energy spectrum is interchanged. 
For the present parameter set, the critical line is shown 
in Fig. 2. When initial couplings {K, A) are located near 
the critical line, the renormalization flow passes near the 
NFL unstable fixed point, and then arrives at the Fermi 
liquid (FL) fixed point. Namely, the NFL behavior can 
be seen in the observable temperature range, if (K, A) 
are close enough to the critical line. Although we present 
here the results for only one parameter set, we have con- 
firmed that the NFL behavior can be seen rather ro- 
bustly at finite temperature even if the conditions, such 
as particle-hole symmetry, spin rotational symmetry, and 
identical Tk's in each orbital, are not satisfied. 

Another fundamental assumption of the present the- 
ory is that Kondo temperatures of two channels are well 
separated as Tki ^ Tk2 due to the difference of hy- 
bridization (Tki = 1.50 X 10-1 and Tk2 = 8.27 x lO"'' for 
J± — Jz — 0). Thus, below Tki, has been screened 
out by conduction electrons in channel 1, while S'2 is 
still active as the localized moment. As the temper- 
ature decreases further, the system sets in the critical 
region, where various quantities exhibit the NFL behav- 
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iors. The characteristic energy Tk is given by that of 
the S=l/2 two-channel Kondo rapdcLiptiuithe exchange 
coupling J - \V2\^/max{Ef2, [/aPBElEf because the 
NFL behavior arises from the balance of two couplings 
by whichjiSgJiiteracts with the two "conduction" electron 
channels £3^ Below the crossover temperature Tcr at 
which X saturates and t^^{uj) ceases to change simulta- 
neously (as shown later in Fig. 4), the system flows into 
the local singlet ground state. Namely, Tcr measures the 
"distance" from the critical point. In fact, as the param- 
eter set {K, A) goes away from the critical boundary, Tcr 
increases. We have checked this by series of numerical 
calculations. Note that the parameter set locating op- 
posite side of the critical boundary leads to the Kondo 
singlet ground state. In the present case, it is estimated 
that Tk ^ 10^'^ and Tcr ~ 10^^, which is assumed to be 
the case of Thi_3;U2;Ru2Si2. The renormalization flow in 
terms of the phase shift is schematically shown in Fig. 3. 

We show the relation between the magnetic suscep- 
tibility, x{T), and the total scattering rate at T = 0, 
r-i(w), for K = 0.14 and A = 0.08 in Fig.4. The tem- 
perature dependence of the resistivity may infer from the 
T~^{uj), which is calculated by sum of each channel con- 
tributions = |I4nP^m('^), being the single- 
particle spectral function. The effect of magnetic field 
has been taken into account through the Zeeman cou- 
pling — {gfiB)jzHz of each /-orbital with j — 5/2 and 
g = 6/7. For Tk < {T,uj) < T^i, x(T) shows the 
Curie- law behavior, while t^^{ll)) takes a broad maxi- 
mum of the order of the unitarity limit. Below Tk, x(T) 
increases logarithmically and t^^{uj) decreases. The 
crossover to the local singlet ground state has not been 
observed in Thi_a:Ua;Ru2Si2, which suggests that Tcr is 
well below the lowest temperature. Behaviors shown in 
Fig.4 well simulates those observed in Ri_a;Ua;Ru2Si2 
(R=Th, Y and La).D^ It is emphasized that the adopted 
K = 0.14 and A = 0.08 is not located so close to the 
critical line. 

The contributions of each channels t^^^{uj) to t^^{uj) 
is shown in Fig. 5. The rather flat bump of t^^{ll)) in 
the region lu < Tki arises from the cancellation of de- 
creasing T^^{uj) and increasing T2^{uj). The gradual de- 
crease of Tj"^ can be understood by an exchange field 
brought about by growing renormalized Hund's-rule cou- 
pling, which destroys the Kondo singlet in channel 1. 
The further flat bottom of t^^{uj) below Tcr indicates 
that the system flows to the local singlet ground state. 

The temperature dependence of the Sommerfeld coef- 
ficient 7=C/T and the entropy under the magnetic field 
is shown in Fig. 6. By applying the magnetic field of 
lO^^D ~ several Tesla, 7 is suppressed at low tempera- 
tures T < Tk because 7 shows stronger divergence than 
— InT in the case oi H = 0. This divergence is caused by 
a transient phenomenon from a NFL unstable fixed point 
to the FL fixed point. The change of the entropy due to 
the magnetic field is moderate compared with that of the 
two-channel Kondo model with magnetic doublet. This 
is because the magnetic field has little influence on the 
CEF singlet ground state. Namely, the behaviors of the 
present model are relatively robust against the magnetic 
field. 



In summary, we have shown that the puzzling problem 
of Thi_3;U3;Ru2Si2 {x < 0.07) can be fully understood 
on the basis of a mechanism characterized by Non- Fermi 
liquid like transient phenomenon with /^-singlet ground 
state. A key assumption is that the system correspond- 
ing to Thi_a;Ua;Ru2Si2 is incidentally located near the 
critical points. However, it is not so artificial to assume 
such an incident if one considers the fact that the NFL 
behaviors appear weakly in much more restricted tem- 
perature range for Ri_2;U2;Ru2Si2 (R=Y and La) which 
would have different set of parameters due to the differ- 
ence of conduction band structure. Experimental results 
show that Tcr in the case of R=La is about lOK while 
that in the case of R=Th has not been observed. This in- 
dicates that the parameter set of Lai_j;Ua;Ru2Si2 would 
be located less close to the critical points than that of 
Thi_j,Uj;Ru2Si2. If Tk's are increased with unchanged 
CEF level structure, the final fixed point can be changed 
from the CEF singlet to the Kondo singlet. Experimen- 
tally, when high pressure is applied, these materials may 
show the increase of the resistivity with decreasing tem- 
perature corresponding to the change of the fixed point. 
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Fig. 4. Temperature dependence of the magnetic susceptibility 
x(T) (triangles) and frequency dependence of the scattering rate 
T~^(w) at T = (squares). 



Fig. 1. CEF level scheme in /^-configuration adopted in this pa- 
per. 
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Fig. 2. Phase diagram of the ground state in A'- A plane. Squares 
represent the critical points determined by NRG calculations. 
The cross represents a parameter set used in this paper. 
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Fig. 5. Contribution of each channels t-^ ^{^) (circles) and 
T2^{'jj) (triangles) to the total scattering rate t^^(uj) (squares). 
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Fig. 3. The renormalization flow below T < Tki ■ The NFL be- 
havior should be observed when the system passes near the crit- 
ical point. 



Fig. 6. Temperature dependence of the Sommerfeld coefficient 
7 = C/T under the magnetic field of H = 10~* (squares) and 
H = 10~^ (circles). The inset shows temperature dependence of 
the entropy. The other parameter set is the same as shown in 
Fig. 4 and Fig.5. 



